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CODED LABEL INFORMATION EXTRACTION METHOD 



2 Introduction 

In certain types of human and machine-readable information-bearing label, a set of 
elements is used to represent the information contained in the label. This 
representation may be by varying the characteristics of the elements comprising the 
label, and also by the position in which the label elements are placed. Reading 
apparatus senses the characteristics and placement of the elements in order to 
decode the information contained within the label. The elements in the label are 
sequentially scanned, to discover the presence of the constituent elements and to 
measure their characteristics and position and thence to decode the information 
contained by the label. 

This application describes label decoding methods by which information represented 
on such a label can be decoded. Typical label embodiments are where the label is 
manufactured from a plurality of magnetically active elements supported on a 
substrate. Information is coded by controlling the relative positions. Independent 
control of the physical properties of the elements, such as shape may also be used. 

In alternative embodiments, information may be coded by controlling the relative 
angles of orientation of the easy axis of the elements, in some cases together with 
independent control of element shape, or other physical property. 

3 DESCRIPTION OF THE INVENTION 



3.1 Scanning Method 

In certain types of label decoding system a label consisting of a set of elements is 
scanned, that is the elements are activated in turn by applying a signal to which the 
scanned el emen t r esponds, aud tu wMc h^he-readingapp a i atu s i s s e nsitive. Either 
or both of the means of generating the scanning signal, and the means of generating 
the received signal in the reading apparatus select a single elements in turn so that 
the signal from a particular element may be resolved from that of the other elements. 

TTremeansofscannmgmaybeautomaticbydectronic^ 

signal and the receive sensitivity along the label or by an automatic mechanism 
causing the label elements to be interrogated in turn. Alternatively the scannmgmay 
be manually operated where the label is moved past the reading apparatus by hand, 
or alternatively the reading apparatus or part thereof may be moved past the label 
by hand. 



3.2 Signal received from reader 
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Fie la shows a plan view of a label where the elements are made out of high 
permeability, low coercivity magnetic thin film material such as Atalante 
manufactured by IST(Belgium). In this case the elements are lmmxl.5mm.The easy 
axis of the material is arranged to be parallel with the longitudinal axis of the label. 

The reading antenna is passed by the label subjecting each element in turn to a 
region of null magnetic field formed by a permanent magnet arrangement A 
superimposed high frequency field, in this case 6.348 kHz causes the material 
element to switch when it is present in the magnetic null. A receive cod is used to 
detect the harmonic signals resulting from the element switching, and the received 
harmonic signal is mixed to baseband and filtered prior to digitisation. Typical 
signals received from such a label, detecting the second harmonic of the element 
switching signals, when the read head was passed in proximity to the teg are shown 
in Fig lb The signal from a single element is shown in Fig lc This shows two half 
peaks of opposite sign with a zero crossing. It can be seen that the label signal 
approximates a supeiposition of single element responses at each position in which 
there is an element. 

The signal is sampled with a fixed sample rate, and converted to a digital 
representation in an analogue to digital converter. The sample rate must be fast 
enough to capture the signal information in accordance with the Nyquist criterion, 
preferably with anti-aliasing filtering prior to the sampling operation. 

The sample number of the centre of each element response depends upon the time 
at which the element passed the read head. 

3.3 Velocity estimates from signal 



In addition, the width of the element response depends upon the speed at which the 
elementpassed the read head. Information aboutthe dynamics lof the i^atove motion 
of the tag and the read head can be inferred from measuring the width of the signal. 
- Th is can be achieved by measuring th e width of t he p o si tive goin g and n egati ve go i n g 
peak associated with the element It can also be achieved by measuring the 
separation between the positive and negative going peaks. This speed estimate is 
essential to estimate the dynamics of motion and to decode the reformation in the 
label reliably without compromising the available code space by adding regular 
features in the label. 

In alternative implementations of the reading apparatus, the scanning motion of the 
magnetic null relative to the label may be achieved by electronically scanning using 
a set of excitation coils that are driven in a phased manner so as to cause the 
magneticnuUtopropagatemspace.ThenuUmaybescannedatsuffic.entspeedthat 

theTwitching of die label elements may be detected direcdy from the receive cods 
withoutneedingto superimpose any separate high frequency excitation. In thwease 
the material element response coupled into the receiver is a single peak. Fig Id 
shows 3 adjacent peaks of the fundamental signal from 3 magnetic dements 
recorded from a reading system employing electronic scanning. In this case the width 
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nf the res0 onse from each magnetic element can be used to provide an estimate of 
tevSVfTTcznnins nuU as it passed over the element being scanned. Tins 
Somation can be used to estimate the detailed dynamics and errors of the system 
a^dXfrommisaHgnmentofmelabelfromthenommdreadaxisof 
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3.4 Amplitude Coded Signals 

Information can be encoded into the label by varying the height of the maternal 
elemTnts, that is the size of the element normal to the longitudinal axis of the label, 
or the width, that is the size of the element along the longitudinal axis. 

Rv maintaining the element width constant, variations in element height 
pLom^Ls^a^ 

2 a for the mechanically scanned second harmonic system, and in Fig2bforan 
decSonicaUy scanned system detecting the fundamental signal received from the 
label. 

a „f linear test labels with 3 elements was manufactured with the element 
Posulns*^ 

Sl wieTmm and the height of the centre element varied between 1mm and 
FoTeach label, in Fig 2a, the peak position, the half peak width, the half peak 

ItivrvaSion of these parameters of the central test dement compared I with Ae 
outer dements, which were maintained at constant hdght, were plotted The 
m«S ureo^ttal ekment amphtude (vampr), rdative to the outer demente, shows 
?S2? nerSnear rdationship to central dement hdght. The measured central 
denStXtive position (vposr), width (vwidr) and separation ^^^ f 
do not vary significantly with central dement hdght, over and above the effects of 
expen^eS eSors and imperfections in the manufacture of the test labds. 

« m ,-1 a rlvFie2b shows therdative variation of measured centrd dement amplitude, 

aan „, bright can be used to code additional infonnation urith a 

be observed. A typical design god is to maximise the range a ^cha Jj" £ 
-d-Acomprom^ 

^enTln oiTTd^a^hdght of labd dements reliably certain of 
&es?d^ente may be constrained to be at a known hdght to act as a reference. In 
Ae exaZk Sated in 3a, the 2 end elements of the label are maintained at 
maxtmZhe. ! ght to act as a reference. A design rule may also constrain a certain 



G:\AMENDED\MJA-FN29 - coded label information extraction method l.wpd 



BEST AVAILABLE COPY 



5 number of the inner elements to also be at maximum height The reference elements 

may be at any a priori fixed position, or at a position determined by an a prion rule. 

The half peak width, or half peak separation can be estimated and used to form an 
estimate of the velocity of the scanning, at the time when the element was scanned, 

10 in the case of a reader detecting the second harmonic of the label signal. 

Alternatively, in the case of a reader detecting the fundamental of the label signal, 
an estimate of the peak width can be used to form an estimate of the velocity of the 
scanning, at the time when the element was scanned. These estimates are typically 
K/w est where K is some constant, and w_est is the estimate of width or half peak 

15 separation, as appropriate. The value of K typically depends on parameters of the 

reader, the element shape and the position of the element relative to the read 
antenna. Since K varies with element properties, to facilitate estimation of the 
scanning dynamics, the process to estimate K may be restricted to elements of the 
same type based on the amplitude of the element signals. 

20 3.5 Peak finding algorithm 

The first step in processing a received signal such as shown in Fig lb or Fig 3b is to 
estimate the parameters of the peaks contained within the signal. This is achieved 
with the following steps: 

25 1 The signal is delimited, that is the continuous received signal is broken up into 

a contiguous set of samples containing the received signal from one label, or a 
set of peaks that is likely to contain enough information to decode the label, in 
the case of non-registered repeating codes. 
2 The signals from this delimited period are sorted into order of ascending 

3 0 magnitude. An indexed sort function, for example using the Quicksort method, 

such as indexx() (Numerical Redpies in C, second edition, ISBN 0 521 43108 5) 
is appropriates© that the original location of the point of any rank can be quickly 
found. 

# 3 The p ositive half peaks a re found by w o rking down from the largest signal in 

3 5 magnitude in the sorted data, testing the signal points as follows. If the test signal 

point is a peak, that is its adjacent signals are lower in magnitude it is estimated 
as a positive half peak, and its amplitude is interpolated by interpolation of the 
peak and its neighbouring points. 

4 From each positive half peak found, the signal is scanned forward and backwards 
40 to find the points at which it falls below a defined magnitude relative to the peak, 

this is interpolated to a fraction of a sample period A typical magnitude value is 
0.8 of the peak value found. The half peak width is estimated as the difference in 
position between these points. The half peak position is estimated as the average 
of these two points. The half peak is rejected if it overlaps with any other pre- 
45 existing peak in the list, that is its width about its position interferes the width 

about any other peak position in the list. If it is not rejected, the half peak is 
appended to a list of half peaks. 

5 The peak finding process stops when the magnitude of the next highest point falls 
below an absolute threshold, or when it falls below a defined fraction of the 
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5 largest magnitude signal, for example 0.3 of the largest magnitude signal, or if 

a limit for the number of peaks is reached. 

6 The negative half peaks are found by working up in magnitude from the lowest 
signal in magnitude in the sorted data, testing signal points of increasing 
magnitude. If the test signal point is a negative peak, that is its adjacent signals 

I o are higher in magnitude it is estimated as a negative half peak, and its amplitude 

is interpolated by interpolation of the peak and its neighbouring points. 

7 For each negative half peak, the signal is scanned forward and backwards to find 
the points at which it increases above a defined magnitude relative to the 
negative peak, this is interpolated to a fraction of a sample period. The half peak 

15 width is estimated as the difference in position between these points. The half 

peak position is estimated as the average of these two points. The half peak is 
rejected if it overlaps with any other pre-existing peak in the list that is its width 
about its position interferes the width about any other peak position in the list 

^ If it is not rejected by this test it is then appended to a list of half peaks. 

20 8 The negative peak finding process stops when die magnitude of the next highest 

point falls above an absolute threshold, or when it falls above a defined fraction 
of the largest magnitude signal, or if a limit for the number of peaks is reached. 
9 The resultant list of positive and negative half peaks is then sorted into range 
order. This process will interleave the consecutive half peaks. 

25 10 Merging range consecutive pairs of half peaks in the range sorted half peak list 

forms a list of double peaks. The amplitude of the merged double peak is the 
difference in amplitude of two half peaks being merged. The position of the 
merged double peak is the average of the two half peak positions. The width of 
the merged double peak is the average of two half-peak widths. The separation 

30 of the merged double peak is the difference between the two half peak positions. 

The results of this half peak finding process are illustrated in Fig 4a which shows the 
half peaks found when the signal shown in Fig 3b was processed, and shows the 
output of step 9 above. These peaks are shown with a cross at the located position 
and estimated amplit ude. Fig 4b shows the magnitude of the d oubl e peaks in 

amplitude units. In this cas e the amplitudes are negative smce the nrsi half peak is 

a negative going peak. The difference in peak amplitude between the full height and 
reduced height elements in the label can be clearly seen. Fig 4c shows the double 
peak separations versus double peak position showing the variation across the 
40 scanning of the label. 

In the case of a fundamental signal where there is just one peak per element, such 
as that shown in Fig Id, a simpler process is may be used, missing steps 6 to 10 
above since there is no need to find and merge double peaks. 

45 10.1 Extraction of velocity function 

There are N label elements positions whose centres are located at X 1 .. X,, along the 
longitudinal axis of the label, and whose position we wish to estimate. It is assumed 
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that X, is at position 0.0 and that the final element, X„ is at die known position of 
corresponding to the length of die label L. 

The peak list as described above contains estimates of the sample number, and hence 
times T, at which the reader detected a double peak. For each element i, the time T,. 
is simply the position in samples where the peak was measured multiplied by the 
sample period, which is the reciprocal of the sampling rate. 

By using the estimated double peak width or double peak separation for each of die 
N peaks denoted W, a velocity estimate for the i* peak is VE,=1/W, at time T, . 

The order M velocity function of time is modelled by a polynomial function. M is 
chosen to allow sufficient degrees of freedom to model the dynamics of label motion 
suffidendy accurately. It can range from 1 (constant acceleration model) upwards 



30 



35 



• W(t,M)=a 0 +a,t+a 2 t 2 +a 9 t s +....a M t M 

where {a} are a set of polynomial coefficients to be determined. 

25 The squared error for velocity estimate VE, at time T, is E, 2 =(VEi-VF(T p lVI)) 2 

The total sum of the squared errors is ES(M)=E 1 2 + E*+ .....+ E» 2 which is 

minimised by varying the M + 1 parameters a^a, a M This is achieved by the method 

of least squares to find the optimum parameters in the least squares sense. 

Define the matrix X to contain the T, values raised to a given power. For example 
column 0 row i of X contains T|0 and column M row 4of X contains T4™ 1 . Define the 
vector V to contain the velocity estimates VE, . For example the element in row 4 ot 
V contains VE 4 

Some linear combination of the columns in X will be closest to the vector V in die 
least-squares sense. The coefficients of this linear combination will be the 
coefficien t* o f die b e s t-fining polynomial. If A is the v eUoi o f length M 1 1 that ho l ds 
the coefficients, the problem is to find the value for A that gives die best near- 
solution" to the system of equations XA= V. If the columns of X are independent, the 
best choice for A is given by: 

a = fx'xy^fx'vi 

Where « T » represents matrix transpose and " " represents matrix inversion. The 
coefficient a, is then the element from row i of A. This may be efficiendy calculated 
by using methods such as LU decomposition, for examples the functions ludcmpO 
and lubksb() (Numerical Recipes in C, second edition, ISBN 0 521 43108 5) 

This process is illustrate in Fig 4d where the velocity estimates W, that were 
computed from the double peak separation data of Fig 4c are shown, those points 
marked "full data" contain the velocity estimates from all the peaks. In this example 
case, only the estimates from peaks with Ugh amplitude, see Fig 4b, ^elected for 
velocity fitting. These selected points are marked "Data points used in fit . If the label 
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known a priori to have dements of the same size then all elements would have been 
used. A least squares fit of order 2 was used and the continuous line shows the 
estimated velocity function as a function of sample number. A least squares fit of 
order 1, 2 or 3 is typically sufficient to model the velocity function with sufficient 
accuracy. 

Error metrics are computed to allow the decoding algorithms to assess the quality of 
the read. A useful normalised velocity error metric, denoted evel, to quantify how 
well the velocity function fitted die data observed is defined. It is the sum of squares 
of the normalised error defined as : 

evel = sum ( EjVVE, 2 ) for i=l..N 

The vector A could also be estimated by an iterative minimisation procedure such as 
Powell's method. 

The above procedure can be repeated for varying model order M on the same data 
in order to find the best approximation. Preferably this is the solution with the 
smallest number of degrees of freedom, thatfits the data read by the reader with the 
lowest error metrics. 



10.2 Analytical integration and normalisation of tag length 

The velocity fit function is VF(t,M)=a 0 +a 1 t+a 2 t 2 +a 3 t a +....+a M t M This may be 
analytically integrated to find the positional fit function 

P(t,M) =Q+a 0 t+(l/2)a 1 t 2 +(l/3)a 2 t 3 +....+(l/(M+l))a M t M+, 
where Q is an arbitrary constant of integration, set to zero. 

We know that the position of label element 1 is at 0.0. We also know that deposition 
of element N, the last element scanned is L, the length of the label. The final 
estimates of element position are: 



XE, b L . (P(X,, M)-P(X„M)) / (P(X N ,M)-P(X 1 ,M)) 



40 
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10.3 Binning of data 

There are N positional estimates XE X to XE N from which N-l gaps may be estimated. 

Typically the gaps have the minimum gap subtracted, and are divided by the 
minimum gap, to give areal number representing the number of gap increments. The 
nearest integer is then taken to represent the gap for subsequent decoding. This can 
be expressed as: 

The N-l normalised gaps are calculated by 

R j =(XE j+1 -XE j -Minimum_Gap)/Gap_Increment for J-1...N-1 
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The nonnalised gaps are then converted to their nearest integer representaton 
Gj=2.t^+0.5) forj=l...N-l.whereint(x) means the truncated integer portion 

of x. 

One useful error metric, denoted emax, is defined as the maximum magnitude of the 
difference between G } and Another useful error metric, denoted eav, is the average 
magnitude of the difference between G, and Rj. 

Fig 4e shows the gap sequence data R, computed from the data in Fig 4a-d, and 
IhoTs the close correspondence or "snapping" to an integer gnd, giving a > fagh 
confidence thatthe gaps have been correcdy measured. The gap sequence in this case 
is 3,2,2,1,7,0,0,2,1 . 

10.4 Enumeration algorithms 

Depending on the coding scheme used, the gap sequence is enumerated to a unique 
decoded lalue. Various codes exist that have different properties for exampkv diose 
codes allowing read direction to be determined, or those that allow non-registered 
TpefattrnXt is allowing the gap sequence to start at any position, for example a 
Portion of a repeating gap sequence. The code may be determined a prion, or 
Sosed by exaniinhig die read signals for the number of material elements 
preSror by estimating the sum of the gap numbers. It is preferable to detenmne 

data. 

For certain codes, the number of elements in the label may be one of a range of 
possibles, with fixed label length in order to increase the number of possible 
30 codes. 

The direction finding codes are preferably used where ^P^^^ ™° d t0 
allow the amplitude information to be associated with the correct element 
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10.5 Amplitude threshold estimation of elements 

Where it is desired to detect the size of the received element signal, for example 
whTthe dements may be of two or more sizes, reference elements may be placed 
I Z set of dements used to make the label. In one embodiment, the end elements 
Se constrained to be the largest element size to act as a reference. Other 
rbodTm^Smayplacereference elements diat have some unique characteristic or 
may be identified by an a priori rule. 
The amplitudes are measured and used to compute normalised amphtudes/This is 
b^uJe P S 

beTe to tne scanning process varying the distance or angle of the read head from 
45 the label during the read, or other factors. 
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The estimates of element position XE, may be used for linear ******* of 
reference amplitude along the label. The reference amplitude predicted at the 
pXon of thedement divides the measured element ampBtodes to 
normalise {he amplitudes. The normalised ampUtudes should be nominally 1.0 for 
an element the same size as the reference element 

Threshold levels are set to provide disoimination between the possible element 

aTpSes and to classify the amplitude of the element into one 

values. For example if there are two values then the amphtude is either high or low. 

JoVexample a half height element may have an amplitade, on average half A e 

reference amplitude ie 0.5. In this case, the thres hold level ^ould be set to 

approximately 0.75 in order to determine if the element was full height or half 

height 

Preferably a code is used which enables the direction of lading of the label to be 
Sermtoed by the sequence of gaps that resulted, to *is cas e i me elements c^i be 
uni^ely identified and it is known which element of the label to associate with the 
decoded amphtude information. 

in the case of the label data shown in Fig 4b, denoting the amphtodes ofAe signak 

from thTdements as «H" for high amphtude and "L» for low amphtude, the elements 

a^go^Ta^ 

b^uSgTbdwith™^^ 

converted to decoded data. 

A typical scheme allows the inner elements only to contain dam igno^ the ^d 
Vlllnt* which are references, and maps H to a logical 0 and L to a logical 1 so that 
Saryd^^^^ 

t 256 Amplitude codes can be represented from 8 inner ^^^^ 
codes used are constrained to ensure that a certain number of elements in the label 
Z ■ exaTpleTare at full height to facilitate velocity estimation. If this constramt is 
useoX number of amphtude codes is reduced by 37 to 219 in this case. 



• An amphtude error metric eamp is d ennea-whidris-th^ 

Ifff^ce between the largest half height element diagn osed to the 
40 heieht element diagnosed. If there is no half-height element diagnosed, eamp is 

sSd elements, a large value for this metric suggests a dear discrimination in 
ZpU^eSieen elements of different types. Smaller values indicate increasmg 
risk that an dement may be categorised as the wrong amphtude. 

45 10.6 Error Metrics 

The error metrics defined: evel, emax, eav and eamp can be used to determine the 
JeadTabr ev2 emax and eav are preferably small and eamp high Independent 
toelhdt for nfese quantities are independently established. These thresholds are 
SleLS values for many labels in many read operations over the range 
50 of reading ranges and scanning dynamics expected m an application. This is done to 
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minimise the probability of a wrong code at the same time as maintaining a high 
probability of a correct read. 

During the read decoding process, if the computed error metrics are higher than the 
threshold value in the case of evel, emax and eav then the read is rejected. . If the 
error metric is lower than the threshold in the case of eamp then the read is rejected. 
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Fig la: Linear Label 
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Received Signal from Label 
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Fig lb: Typical received second harmonic signal from label 

Received Signal from Single element 
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Fig 1c: Typical received second harmonic signal from element of label 
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Fig Id: Received fundamental signal from 3 elements in electronically scanned 
reader 
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Relati ve parameters with element height 
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Fig 2a: Variation 
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Fig 2b: Variation 
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height coding 
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Fig 3b: Receive signal from linear label with height coding 
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Fig 4a: Detected half peak positions 
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Fig 4b: Double peak amplitudes versus tune 

Element amplitudes categorised : H,L,H,L,H,L,H,L,L,H 
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Fig 4c: Double peak separation versus tune 
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Fig 4d: Velocity fit to selected points versus time 
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Number of gap increments 

1 0 Fig 4e.* Resultant fit of Gap number to integer grid 

Nearest integers to gap number is 3,2,2, 1^7,0,0,2,1 
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